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Abstract 

We calculate the second virial coefficient of anyons whose wave function does 
not vanish at coincidence points. This kind of anyons appear naturally when 
one generalizes the hard-core boundary condition according to self-adjoint ex¬ 
tension method in quantum mechanics, and also when anyons are treated field 
theoretically by applying renormalization procedure to nonrelativistic Chern- 
Simons field theory. For the anyons which do not satisfy hard-core boundary 
condition, it is argued that the other scale-invariant limit is more relevant 
in high-temperature limit where virial expansion is useful. Furthermore, the 
cusp existing at the bosonic point for hard-core anyons disappears in general; 
instead it is shown that a new cusp is generated at the fermionic point. A 
physical explanation is given. 
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In this letter we would like to calculate the second virial coefficient of a system of anyons 
1]^ obtained by generalizing the usual hard-core boundary condition at coincidence points 
according to the self-adjoint extension method 0], and discuss its interesting physical prop¬ 
erties. To motivate our calculation, let us briefly summarize the present status. 

Anyons are particles exhibiting fractional statistics characterized by a statistical pa¬ 
rameter a which interpolates between bosonic and fermionic cases. Quantum mechanically, 
they can be regarded as flux-charge composites and their interaction is essentially Aharonov- 
Bohm type [§]. From the held theoretic point of view, they can be obtained if the Schrodinger 
held is coupled to a Chern-Simons gauge held @,||. The connection of these two descriptions 
is now well established and an interesting subject here is the relation between the method 
of self-adjoint extension in quantum mechanics and the regularization-renormalization pro¬ 
cedure in held theory where |0|^-type contact interaction is necessarily introduced l^-Q. 
Recently, an all-order analysis has been carried out to establish the precise correspondence 
and the role of contact term has been emphasized in ref. [Q. 

In the literature, regularity of wave functions is often assumed and it leads to the hard¬ 
core boundary condition |^. Although this hard-core case is conceptually simple and easier 
to treat, there is no a priori reason to assume such a boundary condition in general; a real 
system under study should determine the relevant boundary condition eventually. If we 
relax the regularity requirement, allowing wave functions to diverge at coincidence points 
according to the method of self-adjoint extension, we get a one-parameter family of boundary 
conditions which in general introduces a scale in the theory. (The corresponding anyons were 
called “colliding anyons” |^.) The scale invariance of the theory is restored only in the 
limit that the scale parameter goes to either zero or inhnity. Field theoretically, when the 
strength of the |(/)|^-interaction is equal to the critical value A = = ±Ao, the held 

theory becomes ultraviolet hnite and scale invariant [Q. These critical values have some 
more interesting properties. The repulsive case A = Aq corresponds to the usual hard-core 
boundary condition, while in the attractive case A = — Aq, it is known that the classical 
theory admits static soliton solutions satisfying a self-dual equation p. Also, for A ^ — Aq, 
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after performing renormalization, one can relate the renormalized coupling Aren to a nonzero 
finite value of the self-adjoint extension parameter in a specihc way [pd],p|. 

Since the free anyon gas has not been exactly solved, one usually resorts to the virial 
expansion to study thermodynamic properties in high-temperature low-density limit. For 
example, the second virial coefficient was calculated for anyons with hard core in ref. |T^-jI5[ . 
Here an interesting fact is that, according to recent analyses [01, Aren = Ao (hard-core case) 
is an infrared hxed point, while Aren = —Ao is an ultraviolet hxed point, with Aren flowing 
from Ao to —Ao as the renormalization scale increases. This implies that the virial expansion 
for usual anyons with hard core is not much relevant for the case of colliding anyons in high- 
temperature limit where the virial expansion is useful. Therefore it may be interesting to 
calculate the second virial coefficient for anyons with an arbitrary boundary condition and 
study the behavior as the scale changes. Indeed, we will see that the result of hard-core case 
does not represent a typical behavior of the other cases. For example, the cusp existing at 
bosonic point a = 0 in hard-core case is a very special feature of scale-invariant limit which 
does not occur in general; instead a new cusp is generated at the fermionic point |a| = 1 
except for the hard core case. We will give a physical interpretation of this result later. 

Here we should mention earlier works 100 which have a partial overlap with this 
paper. They calculated the second virial coefficient for the case that the system has a 
bound state. However they did not pay much attention to physical properties for general 
boundary conditions except at the hxed point, Aren = —Aq. Also, we believe that our method 
is simpler and physically more transparent. 

Let us start with the Hamiltonian of a system of free anyons with mass m and statistical 
parameter a in bosonic description 




( 1 ) 


where 


Ain = E 

m(^n) 


Xt — Xi 


( 2 ) 
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is the Aharonov-Bohm-type vector potential and r = denotes the position of the 

n-th particle. As mentioned above, this Hamiltonian does not specify the system completely 
and one has to give a suitable boundary conditions when particle positions coincide; we will 
give the form later. 

In two dimensions, the second virial coefficient B{T) is given by 


B(T) = A(--Z^IZl), 


( 3 ) 


where A is the area of the system and Zi and Z 2 are the one-body and two-body partition 
functions, respectively. They are calculated as 


Zi — AX^ , 


Zo = 2AXI.Z, 


T-^rel 1 


( 4 ) 


where \t = {271 /rnksTyB jg thermal wavelength and Z^^i = Tre is the single 
partition function of the relative dynamics with 


^rel = — (P - aA)^ 
m 


( 5 ) 


After the separation of variables V’(r) = e^"'^Rn{r), the Schrodinger equation for H^-ei 
reads 


1 

m 


Id d (n -f- ay 
r dr dr 


^2 

Rn{r) = ERn{r) = —Rn{r) , 

m 


( 6 ) 


where n is even since we are working in bosonic description. Also, since a enters only 
through the combination n + a, we may restrict a to the interval [—1,1] without loss of 
generality in considering the energy spectrum |l|. Equation is nothing but the Bessel 
equation and the general solution is of the form 


Rn{y (/cr) (/cr), 


( 7 ) 


where A, B are constants. For n 7^ 0 the constant B must be chosen to be zero for square 
integrability, but for n = 0 (s-wave case) arbitrary A and B are allowed. This leads to an 


above-mentioned one-parameter family of boundary conditions for s-wave | 19 ], 
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( 8 ) 


with the corresponding solntion given by 


i?o(^) = (const.) 


2|«l 


+<^ I -| .^-|n| (*:’■) 


(9) 


where a = ±1 and k is a scale introduced by the boundary condition. If k —> cx) with 
a = +1 we recover the hard-core boundary condition '0(0) = O.Q In addition to the solution 
(I), if cr = — 1, there is a bound state satisfying (H), 


/?o(r) = (const.)i^|„|(Kr) , 


( 10 ) 


with energy = —K^/m. (Thus the parameter k in the boundary condition appears in the 
bound state energy.) To proceed, we introduce a circular boundary at radius R and assume 
Ro{R) = 0 as usual 0. The raduis R will be taken to be infinity eventually. Then the 
allowed energies are 

kl 


k^n,s 


(n ^ 0), 


m 

Eo,s = - (n = 0), 


m 


where kn^sR is the s-th zero of J\n+a\{kR) = 0 and kgR is the s-th zero of (^, i.e., 

/^\2|“l 

J\a\{kR) + a i-j J-\a\{kR) = 0. 

Hence, compared with the hard-core case, only s-wave energy spectrum is changed. 
From Eqs. (^) and (^, the second virial coefficient can then be written as 


B{T) = Bq{T) - 2X1 { e-^^^e{-a) + lim V 

I s=0 


3 -/ 3 - 


3 -/ 3 ' 


Q,s 


( 11 ) 


( 12 ) 


(13) 


^If a = —1, formally we obtain 0(0) =0 when k —> oo but it does not correspond to the usual 
hard-core limit. Actually, the system is not bounded from below because it has a bound state with 


negative infinite energy as seen in (IT) below 
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where (3 = l/ksT and Bq{T) is the second virial coefficient for hard-core anyons calcnlated 
in ref. |T^ , 


Bo{T) = A| 



(14) 


To perform the snmmation in (|^ in the limit R 
cos7rQ;J|Q|(a;) — sin7r|Q!|A^|c|(x), 


cxo, we rewrite (|^ nsing J_|„|(x) = 


cosri{k)J\a\{kR) — sinri{k)N\a\{kR) = 0 , 


(16) 


where r]{k) is dehned as 


taiiT]{k) = 


/fc\ 2 |o| . I I 

CT1-1 sm7r|Q!| 


,2|a| 


COS 7ia 


i+ht) 

Then, from the asymptotic behavior of Bessel fnnctions 


(16) 





X — [\a\ + 


1\ vr 


2 / 2 


1 \ TT 

^-ll«l + 5) 5 


+ ■ • • ) 
+ ■ ■ • ) 


(17) 


it is easy to see that, for any e > 0, there exists an integer sq independent of R snch that 

< sin e , if s > sq- (18) 


cos 


ksR- ( |a| + ^ + r]{ks) 


In other words, if s > Sq, 

ksR + ri{ks) = STT + +0{e). (19) 

Now we split the snmmation into two parts Y3T=o = • The hrst term will 

vanish eventnally ns R —>■ oo and so the second term dominates the whole snmmation. For 
s > So, we may approximate fcg’s in the snmmation by (|T^ (omitting 0(e) term) and the 
error of the approximate snm is easily seen to be 0(e ■ R^), which vanishes as e —0 and 
R ^ oo. This allows ns to replace kg by a fnnction k{s) dehned for real s as 

k{s)R + r]{k{s)) = S71+ (20) 
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when s > sq. [k{s) is defined similarly.] Then by use of Euler-Maclaurin formula ||^, the 
summation may be converted into an integral, 


E 

s=0 


/3fcf 


/3fc, 


O.s 


C ^ ^ 


ds 


f So 


+ 0(e-R’')+0{l/R), 


( 21 ) 


where we did not write explicitly other terms in the Euler-Maclaurin formula which are 
0(1/R) at most. From (pOD, ^ = ^{R + and % = Thus we get0 


lim 'Y' 


R-^oq 


s=0 


dk 

C ^ ^ 


1 dr] 1 

- / dk-^e -f- - 77 ( 0 ) 
TT Jo dk TT 


( 22 ) 


From (^), 77 ( 0 ) = 0 unless k = 0 and so we drop this term from now on assuming k 7 ^ 00. 
^ can easily be computed from ([Tbl) . Inserting this into we finally obtain 


B{T) = Bo{T) - 2X1 \ e^e{-a) + sin 7 r|a| / 

TT JO 


0 1 + 2a COS 71 at\^\ + 


(23) 


where e = fdK^/m. If we put a = —1, this result is equivalent to the formulas in refs. ||T6|,p!7 


The above integration cannot be performed analytically in general. In case of semions, 
|a| = 1 / 2 , it is reduced to the error integral. 


B(|a| = 1 / 2 ,T) = B(,(|a| = 1 / 2 , T) - A|e'|l - cTErf(v^)], 


(24) 


where Erf (x) = fo dt. Compared with i?o(|tt| = 1/2, T) = A^/S, the second term 
becomes dominant for most of £ except near the hard-core limit e ^ oo with a = +1. 

For general value of a, one can obtain various limiting behaviors by examining the 
integral. First, for e: S> 1, we find the asymptotic behavior 


B{T) = Bo{T)-2Xl 


e^9{-a) 


crr(|a|-hl) . 


TT 


r|a| 


SlUTT a -t- 


(26) 


We see that, for cr = -|-1, the hard-core result is recovered but, if a = — 1, the bound state 
contribution dominates the virial coefficient. Because of the nature of the virial expansion. 


^With the second term Ar/(0), (22) generalizes the usual formula in terms of phase shifts |21| 
whose naive application does not produce the correct result in anyon case as noted in 
^When AC = 0, 7y = 7r|a| and the first term involving ^ becomes zero. 


7 































(p5|) is physically meaningful only for parameters in the region K^/m S> ksT p/m where 
p is the particle number density. For e 1, which corresponds to high-temperature limit, 
we hnd after some calculation 


B{T) = Bo{T) - 2A||a|(l - ae'"' + ■■■), 

= Ay - — |a| — -|- 2(TAy|a|£l"l + • • • , (26) 

which differs from the hard-core result by —2Ay|a| for e: = 0. 

One may also consider the behavior of B{T) near the bosonic {a = 0) and the fermionic 
(|a| = 1) points. After some calculation, we obtain^, for a 0 and for nonzero hnite e, 


B{T) = Bo{T) - 2X^u{e)9{-a) - |a|A^ + 0{a^) 


2u{e)9{-a) 


Ay -h O(a^) , 


(27) 


where is the Neumann function dehned by 

, , /■°° e^dt 


i^[£] = 


(28) 


/o T(t -|- 1) 

Therefore the linear terms in |a| are cancelled out, i.e., B{T) is a smooth function of a near 
a = 0 and the cusp at the bosonic point existing for hard-core anyons disappears. On the 
other hand, near the fermionic point, we hnd 


B{T) 


Bo{T)-2\^^e-^^ + Ue){l 


A^ 


■o(l-|«l)' + 7-2e- 


a|) H- 

+ /o-(^)(l 


a|) -t- • • • , 


(29) 


where {o' = ±1) are some complicated functions of £ involving exponential integrals 

whose explicit forms are not much illuminating. Thus we see that in general a cusp exists 
at the fermionic point |a| = 1 for anyons without hard core. 


“^When a = —1 and a = 0 this result gives the second virial coefficient of bosons with the self- 
adjoint extension discussed in ||6|,pT|. 









The physical reason for this change of the cusp position may be the following. For the 
hard-core boundary condition, anyons behave more like fermions in the sense that they 
cannot overlap and the lack of hard core for bosons is responsible for the cusp |^. In 
contrast, for other boundary conditions, anyons can overlap like bosons and this generates 
a cusp at the fermionic point where particles cannot overlap. 

We plot these results in Figure 1. The deviation of the virial coefficient from hard-core 
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result is quite signihcant for reasonable values oie = (5^. Also, Figure 1 clearly shows cusps 
at |a| = 1 and smoothness of curves at a = 0. 

One may be able to generalize this calculation to a system of multi-species of anyons obey¬ 
ing mutual statistics [^. This system is known to be related with multi-layered quantum 
Hall effect. In this case self-adjoint extension is allowed in p-wave part as well as in s-wave 
part and it introduces another scale parameter. Field theoretically it would correspond to 
introducing a nonrenormalizable contact interaction with derivatives of delta functions. It 
may be interesting to see the effect of these scale parameters on various physical properties 
of the system. Also, the analysis of this paper may be able to be generalized to a system 
of non-Abelian Chern-Simons particles for which the second virial coefficient was recently 
calculated in the hard-core case 0- 
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continuous encouragement, and S. J. Kim for useful discussions. We also thank A. Moroz 
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FIGURES 



a 


FIG. 1. The second virial coefficient of anyons without hard core as functions of the statistical 
parameter a for various values of e = j3K^/m. e = oo with a = 1 corresponds to hard-core limit 
and e = 0, high-temperature limit. 
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